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Abstract 


The general formulation of a technically advantageous method to find the 
ground state solution of the Schrodinger equation in configuration space for 
systems with a number of particles A greater than 4 is presented. The wave 
function is expanded in pair correlated hyperspherical harmonics beyond the 
lowest order approximation and then calculated in the Faddeev approach. 
A recent efficient recursive method to construct antisymmetric A-particle 
hyperspherical harmonics is used. The accuracy is tested for the bound state 
energies of nuclei with A = 6 -i- 12. The high quality of the obtained results 
becomes evident from a comparison with other approaches. 
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I. INTRODUCTION 


Few-body nuclei with a number of nucleons A between 5 and 16 are a particularly inter¬ 
esting testground for nuclear theory. They lie in the range between the classical few-nucleon 
systems {A < 4) and the smallest nuclei, which can be described realistically starting from 
a mean field ansatz. Therefore one hopes that these nuclei could build a link between few- 
body and many-body physics. Presently quite an effort is made for a better understanding 
of these intermediate systems. Specific interest is devoted to halo-nuclei, but also the less 
exotic nuclei of this mass range are investigated thoroughly. Many theoretical techniques 
for the calculation of their ground states have been imported from the classical few-body 
held, where there has been a considerable progress in the last years. In fact for the classical 
few-body systems rather different approaches have been developed and proven to lead to 
precise results. These methods include solutions of the Faddeev-Yakubovsky equation, vari¬ 
ational (VMC) and Green Function Monte Carlo (GFMC), the Hyperspherical Harmonic 
(HH) ansatz, the stochastic variational (SVM), as well as coupled cluster and resonating 
group methods. 

For nuclei with A > 4 a similar level of precision has not yet been reached. Here 
GFMC has led to the most accurate results. Exact bound state energies with realistic NN 
interactions have been calculated for A < 9 [Q. Unfortunately wave functions cannot 

be generated with this method. Recently, a very powerful tool to calculate few-body wave 
functions has been developed with the SVM ansatz However, this approach is probably 

most suitable for systems with A < 8. For nuclei with A > 8 rather good results have 
been obtained in the Integro-Differential Equation Approach (IDEA) 0], which uses the 
HH expansion at hrst order, and with a variational method, the Translational Invariant 
Conhguration Interaction Method (TICI) P,|^, which is inspired by the coupled cluster 
method. 

In this work we present the general formulation of the method, which combining the 
main ideas of the HH expansion, the pair correlation ansatz and the Faddeev approach 
allows to calculate wave functions of few-body systems. We apply it to calculate the binding 
energies of nuclear systems with 6, 8 and 12 particles. Including higher order HH functions 
we make the nontrivial step beyond the IDEA approach. The difficulty in constructing 
antisymmetric A-particle HH is overcome by the use of a recently developed very efficient 
recursive method 0, where HH basis functions are constructed, belonging to well dehned 
irreducible representations of the orthogonal and symmetric groups. 

The paper is organized as follows. The method is described in Section |^, while the 
construction of the HH basis function is briefly reviewed in Section |IH[ Section illustrates 
how the matrix elements are calculated and the obtained results for the binding energies of 
nuclei with A = 6, 8 and 12 are discussed in Section 0. Gonclusions are drawn in Section 

B 


II. GENERAL FORMULATION OF THE METHOD 

Our aim is the solution of the Schrodinger equation for a system of A particles interacting 
via a two-body potential. After subtraction of the center of mass Hamiltonian the problem 
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reduces to the solution of the internal Hamiltonian in terms of = (H—1) Jacobi coordinates 
fji. One can write the wave function in terms of an HH expansion 

^{f]i..f]A-l,Si..SA,ti..tA) = '^RKu{p)HKu{^,Si..SA,ti..tA) 

Kv 

= ^^Kv{j]i-,Sx-SAM-iA) ■, ( 1 ) 

Kv 

where, in case of nucleons, 

Si-.Sa, = X! y KNLNMNFAYA-ia^i^) ^ 

ya -1 virAl 

'^SS"^Tr^fyl,’h4_iaST(si..SA, tl-.tu) (2) 

represent the hyperspherical harmonic functions coupled to the spin-isospin basis functions 
to yield the totally antisymmetric wave function. The yKNiNUffTAYA-ia^i^) 
called symmetrized hyperspherical functions that depend on the N Jacobi coordinates. They 
are HH functions with hyperspherical angular momentum K = and good angular 

momentum quantum numbers Ln, Mn that belong to well dehned irreducible representations 
(irreps) Ti G r 2 ... G Ta, of the permutation group-subgroup chain C ^2 .. . C 5a, 
denoted by the Yamanouchi symbol [rA,YA_i] = [Ta, Ta-i, ..., hi]. The dimension of 
the irrep r„ is denoted by |r„| and Kya,Ya-i ^ phase factor [§]. Similarly, the function 
^ss^tt^VaYa-ici^/ ^ symmetrized spin-isospin state with good quantum numbers 5, 5^,T 
and T^. The label (oa^) is needed to remove the degeneracy of the hyperspherical (spin- 
isospin) states with a given symmetry. The argument of the hyperradial function Rkv is 
the hyperradius p dehned by 


P 


= Y.n. 


2 

i 5 


( 3 ) 


H is the (3A-4)-dimensional hyperangle, and for brevity we shall use v for all the quantum 
numbers but K, i.e. z/ = SS^TT^ Tacc^cx^). In what follows we shall use the 

subscript N to denote A particle quantum numbers that depend on the Jacobi coordinates, 
and the subscript A — 1 for the A — 1 particle subsystem. 

It is well known that hnding a solution of the Schrodinger equation in terms of the ansatz 
(0) can be very difficult because the number of basis functions increases very fast with K 
and in order to have a real convergence one must use a huge number of basis functions 0 . 
Therefore a correlation function is advantageous to give the wave function a proper behavior 
0. Its advantages have been extensively verihed for classical few-body systems reaching a 
high level of accuracy [pT| , p!^ . A general ansatz within the two-body correlation scheme is 
the Jastrow factor 


=n /«*. (4) 

i<j 

where fij is a two-body correlation function. However, the use of the Jastrow ansatz leads 
to 3A — 3 dimensional integrals. Therefore it is more convenient to use the so called pair 
correlation ansatz 
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( 5 ) 


I” = E X'.,*. 

i<j 

because in this case one can use the Faddeev approach which leads to at most four-body 
integrals. In the Faddeev approach the Schrodinger equation is replaced by equivalent 
equations, 


{T - , (6) 

where d' = J2i<j In order to speed up the convergence these equations can be further 
modihed 0] to include the contribution of the hypercentral potential explicitly. The hy¬ 
percentral potential I4c(p) is dehned as the projection of the two-body interaction on the 
subspace of the lowest order hyperspherical state, i.e. the hyperspherical state expressed in 
Eq. (^, with the minimal K: 

Vhcip) = j . (7) 


Here the integration is carried over the ?>N — 1 dimensional hypersphere, and an implicit 
summation over all the spin-isospin states is understood. With the help of this dehnition 
we can rewrite Eq. (|^) as follows 




E 


= -[^P - Vhcip)]^ 


( 8 ) 


Motivated by the pair-correlation ansatz we shall expand the Faddeev amplitude in 
the following way. 


mp 'y I^)Xp('^p) 'y ' Rkvp^P^HS\..SA-, tl-'tA)yip{^^ij') 1 ( 9 ) 

P Kvp 

where Xp(%) is a polynomial of order p, and % is related to the relative two-body distance 
through 


Zij = Tij/p. 


( 10 ) 


Although expanding the correlation function in terms of Zij should be equivalent to an 
expansion in it is more stable numerically to take Xp{zij) since Vij goes to zero with p. 
Substituting the expansion of Eq. (H) into Eq. (j^) we get 


1 . A(A - 1) 


^n + 


VUp) - E 


71=1 


II %Xpizij) = -[Vij - Vhcip)] Xpizki) , ( 11 ) 

P P k>l 


where A„ is the Laplace operator associated with the n’th Jacobi coordinate fjn- Multiplying 
Eq. o on the left by iHx'u'Xp'izij))^ and using the notation 




< K'u'p'\d\Kup >kl= J dn Hl^,^,in, Sn,tn)x*p'iZij)OHKuiV., Sm,tm)XpiZkl) (12) 


we get 
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( 13 ) 


1 A(A — I'l 

ij < K'v'p'l - -^- -Vhcip) - E\K up >ij Rxupip) = 

Kvp n=l 


- Y1 ^J<K'u'p'\[Vij-Vhc{p)]J2\^"’^”P>klRK"u"p"{p). 

K"v"p" k>l 


If the norm matrix 


Ak'u'p' ,Kup — ij^ E up \Kup ^ ij ( 14 ) 

is non singnlar, one can define the projection operator 

Q = H \K'u'p > < Kup\ . (15) 

Kup^K'u'p' 

The projection operator commntes with Vij and one can rewrite Eq. ([TB| ) in the following 
form 


1 A(A — 1) 

ij < K'u'p'l - - -z- -Vhcip) - E\Kup >ij Rxupip) = 

Kvp ^ n=l ^ 


-E E E «<AVp'||l/„.-l4o(p)l|A"VV'>«x 

Kvp K"v"p" K"'v"'p"' 

^K"v"p" ,K'"v'"p"' ij ^ K"'u"'p"'\Kup >ki RkvpIp) . 
kl 


( 16 ) 


Using the hyperspherical coordinates, and the following definitions 


Tk^u'p'^K ipp ij ^ ^ ly p \K\Kiyp ^ij , 

Vk;v.K^{p) = « < AVVm^ - l4.|AVp >„ , 
lUL'ptiCi/p = y yK^v'p',K"v"p" ij < R"y'p"\R^P >ik ^ 

K"v"p" 

^\<:'v'p',Kvp = y yK'v'p',K"v"p" ij < K U p \Kup >kl , 

K"v"p" 


Eq. (P!BD can be expressed as an hyperradial equation 




dP 3A-4: d 1 A{A - 1) 


p dp 2 


+ 


-vup)y 


k + i,J , 

k,l ^ iJ, 


pp' 


E A^^i 


Rp'{p) 


(17) 


= - i: q,,(p) (i„>» + 21A - 2)»--,li„ + 2)(^ ^M v') V(p) ■ (18) 

p'p" \ ^ / 


Here p = {Kup), and K is the generalized, hyperspherical, angular momentum operator. 
This hyperradial equation can be solved by expanding the hyperradial function Rp{p) into 
basis functions (f>np{p) 


RM = i:c;y„,(p) 

rip 


( 19 ) 


containing the Laguerre polynomials L“ 
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( 20 ) 


'i’nJp) = 


n 


p- 


-3(A-1) 


\ (rip + a)! 


a-(3A-4) 
2 




exp 


P_ 

'2b 


Multiplying the right and left hand sides by 4>np{p) performing the hyperradial integra¬ 
tion one remains with the generalized eigenvalue problem 


- EMmn) = 0 , ( 21 ) 

m 

where m and n stand for the sets K,i',p,np and K', u', p', n'p, respectively, and the matrix 
elements and M^n are given by the corresponding hyperradial integrals. 

Eq. ([TsD represents the main equation of the present method. The following remarks 
have to be stressed in order here: 


1. the 3(A-l)-dimensional Schrodinger equation has been reduced to a 3 x 4 dimen¬ 
sional integro-differential equation, due to the presence of at most four-body terms 

2. to the extent the convergence in the expansions is reached the result for 4/ = J2i<j 
can be considered as an ’’exact” ground state solution of the many-body Schrodinger 
equation. 


III. CONSTRUCTION OF THE BASIS FUNCTIONS 

Technical difficulties in calculating the matrix elements and are encountered 
in the construction of the basis functions Hku{^i Si^ti) for increasing numbers of particles 
and hyperspherical angular momentum K, as well as in the calculation of the two-, three- 
and four-body terms that appear in and M^n- The basis functions Sj, f*) must 

be totally antisymmetric functions in total space (coordinate, spin, isospin). Therefore, one 
needs efficient algorithms to construct convenient expressions for the symmetrized spin- 
isospin basis functions as well as symmetrized hyperspherical wave functions. 


A. Spin—isospin states with arbitrary permntational symmetry 

The construction of A-particle symmetrized spin-isospin functions as well as hyper¬ 
spherical functions can be done recursively ^j. Each A-particle spin-isospin function of a 
well defined irrep of the symmetric group T^ and spin-isospin quantum numbers Sa,Ta, 
is written as a linear combination of spin-isospin coupled products of an (A-l)-particle 
wave function of well defined irrep r^_i and spin-isospin quantum numbers Sa-i,Ta-i, 
and the A’th particle wave function with spin-isospin quantum numbers SA,tA- The co¬ 
efficients of this linear combination are a sort of coefficients of fractional parentage. The 
subspace of functions with good quantum numbers Sa, Ta and T^-i is an invariant subspace 
for the transposition class-sum operator of the symmetric group Sa- The eigenvalues that 
are obtained after the diagonalization of the class-sum operator identify the irreps of the 
symmetric group uniquely. The eigenvectors are the spin-isospin coefficients of fractional 
parentage (stcfps) 
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isiti) = Y. 

5yi_iTA-ia3i-i 


[{SA-i;sA)SA{TA-i;tA)TATA-iaYi\}SATATAa^A ] x 

[^SA-lS^A-lTA-lT-A-irA-iyA-20^l^^ ® 


( 22 ) 


The sum over Sa-i and Ta-i in Eq. (|^) is subject to the usual angular momentum coupling 
rules. 


B. Symmetrized hyperspherical functions 


In an analogous way, one can use the recursive methods developed in the last few years 
^,n] for constructing hyperspherical functions that belong to well defined irreps of the 
symmetric group. In this methods the reduction problem O 3 A -3 D O 3 ^ Sa is solved and 
one obtains hyperspherical functions which belong to irreps of the symmetric group Sa and 
have good angular momentum and hyperspherical angular momentum quantum numbers. 
These functions are expressed in terms of 5 a-i hyperspherical states, that are coupled, 
via the “tree” method [^, with the appropriate single-Jacobi coordinate hyperspherical 
functions, into coupled A-particle states with the desired angular momentum and hyper¬ 
spherical angular momentum quantum numbers. Actually, these coupled states yield an 
invariant subspace with respect to Sa- The transposition class-sum of the symmetric group 
(the second Casimir operator) is diagonalized within this subspace. The eigenvalues that 
are obtained after the diagonalization identify the irreps of the symmetric group uniquely, 
and the eigenvectors are the hyperspherical coefficients of fractional parentage (hscfps). 

Constructing basis functions in such a recursive way makes the evaluation of any two 
body operator easy, since only the matrix element of the two-body operator between the 
last two particles needs to be calculated. A further improvement in the efficiency of the 
algorithm can be reached if one uses reversed order Jacobi coordinates for the construction 
of the hyperspherical functions. Normally Jacobi coordinates are dehned so that the hrst 
one is the relative distance between particles 1 and 2 , while the last is the distance between 
the A’th particle and the center of mass of the A -1 particle system. This implies that in 
order to calculate the matrix element of a two-body operator, depending typically on the 
interparticle distance between the A’th and (A-l)’th particles (see Section |^), one needs to 
rotate the last two Jacobi coordinates by a proper angle so that one of them represents the 
interparticle distance. This can be done by using the Raynal Revai, the T-coefficients and 
the 6 -j coefficients. Constructing the Jacobi coordinates in reverse order, would simplify the 
calculation of two- and three-body matrix elements as no further rotation is needed [T^ . 


Recently an alternative way to construct the symmetrized hyperspherical functions 
has been proposed |0. This method consists in introducing as an intermediate sub¬ 
group the orthogonal group of kinematic rotations Oa-i, i.e. one uses the group chain 
O 3 A -3 D O 3 0 Oa -1 D O 3 0 Sa- The introduction of the kinematic group turns out to 
be of particular importance for increasing number of particles and for large values of K as 
it results in a signihcant reduction in the number of hscfps and in the computation time. 
Another benefit in using this method is the realization of kinematic rotations through the 
representation matrices of the group Oa-i thus avoiding the use of the Raynal-Revai and 
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the T-coefficients in calculating matrix elements of two-body operators depending on the 
interparticle distance. In the present calculation we have adopted this procedure, even if 
it is less efficient than dehning Jacobi coordinates in reversed order. The reason is that in 
this way one can easily extend the present method to solve the Schrodinger-like equation 
with source, necessary for the calculation of the Lorentz integral transforms of response 
functions [0. The presence of the source annuls the advantages of the alternative set of 


Jacobi coordinates. 

Summarizing the procedure for the construction of the symmetrized hyperspherical har¬ 
monics one can split it into two steps. In the hrst step one constructs hyperspherical func¬ 
tions with good orthogonal symmetry, i.e. basis functions with good quantum numbers 
Kn, Ln, Mn that belong to a well dehned Gel’fand-Zetlin pattern = (Atv, An-i, ■ ■ ■, A 2 ). 
The orthogonal group irreps Xj,{j = 2, 3,...) are characterized by the integer or half integer 
numbers Aj^i, \j^ 2 , ■ ■ ■ , ^j,k, where k = [|]. In Ref. [Q it was pointed out that for hyper¬ 
spherical functions these numbers are always integers and, in addition, there are at most 
three non-zero values in the irrep Xj for j > Q i.e., Aj,i, Xj^ 2 , and Aj^s. The second step is 
to reduce each irrep of the orthogonal group On into irreps of the symmetry group Sn+i- 
These two steps are carried out using the recursive method developed in Ref. 0. At the end, 
the symmetrized A-particle, A^-Jacobi coordinate, hyperspherical states can be expressed 
in terms of (A — l)-particle states coupled to the A’th particle state by using two new types 
of coefficients of fractional parentage, namely the orthogonal-hyperspherical cfps (ohscfps) 
for the construction of hyperspherical functions with good orthogonal symmetry and the 
orthogonal cfps (ocfps) for the reduction Oa-i i One has 

^iCjvnivMjvAAr/3|^rAYA-i 

= X! [(AAf-ir^-l/^A-l) X 

(-^N-l-kN-lXN-l/^N-lJ ^N)-^N-kNj}-^N-kNX 

y(Kjv-iLjv-iXjv-i/3^_,rA-iYA-2/3^_,-/jv)KjvLjvMjv ' 


K 


(23) 


Note that the order of the summation in this equation is important. We have to start from 
the sum over the ocfps 

[(A^,_ir^_i/?i_i) A^|}Ajvr^/3i], (24) 

to determine the irrep A^v-i of On-i and then we can sum over the ohscfps 


{Kn-iLn-iXn-i(3n-i', iN)KNLN\}KNLNXNf3N 


( 26 ) 


where {j3§ and /3^) are the degeneracy removing labels {{Pn^Pa) = ^n)- 


IV. CALCULATION OF THE MATRIX ELEMENTS 


The matrix elements of any two-body operator 









Oij = 0^{fij)0^'^{si, Sj, u, tj) (26) 

between the fermionic, antisymmetric, hyperspherical basis functions of Eq. (j^), can be 
written as a sum of a spatial term multiplied by a spin-isospin term. 


< Ku\Oij\K'u' >= E n( ISF 


X 


rA-irt_,ryi_2 


< SAS^ATAT^ATAYA-ia^A^\0^^{s,,s,,U,t,)\S'^S'^AT'^r^AT'AYA-iCyT > X 


where the sum runs over all the 5^-2,5^-1 irreps rA-ir( 4 _^rA -2 such that Ta -2 G E^-i G E^ 
and TA -2 G r( 4 _^ G F'^. The factor preceding the spatial and spin-isospin matrix element 
is a product of the inner product symmetric group isoscalar factors for the antisymmetric 
representation and is given by 


n( ISF ) = A; 


r^TA-l ' 


IFa-i 

- A T2 T2 

|rx-2| , 

IF' 1 
d A-ll A 

(PI 

1 

l-i 

^ ITaI 

A-li A-2 A 


|PJ 

lr'24-ll 


(28) 


The phase factor Ar^PA-i [|I3 is positive (negative) when the number of boxes in T^ below 
the row of the A’th particle is even (odd). 

Using the permutation symmetry of the HH states the matrix elements of any two-body 
operator, Oij, are equal to the matrix elements of Oa,a-i- The spatial matrix elements 
are then calculated using the last two generations of cfps, referring to the constructions 
A — 2 —*■ A — 1 and A — 1 —> A, and the appropriate kinematical rotations which reduce 
the calculation of any two-body operator to sum over one dimensional integrals. The formal 
derivation of the calculation of two-body operator matrix elements between the symmetrized 
hyperspherical harmonics of Eq. (^) has been presented in rather a detailed manner by 
Barnea and Novoselsky |^. Summarizing their results one has 

'^rA-2,r'^_2 E ^ E E 

Ajv-i/3^_]^ Ajv_i/3^7i Ajv-2/3^_2 

A„|}A„r^;3^]*|(A;v-ir:,_i^iy,) AlvDAWtl 

x[(AA,_jr.t-2/3y2) AA,_,|}Ajv_irx-i/3y,]-|(AA,_2r2,-2/3y2)Ay,|}Ayir'^_i/3y,] 

A'„ 


X i: 

\" 


D 


\n * 


i9A,A-l) D 


{.gA,A-i) 


X < iF^LArM^AAr/?^A'^_iU4-2/dE2)|0^(y26v)|i^)vT)vTf;A)v/?^'A';_iU4-2/3E2) > 

Here, for simplicity, we use the symbol A^r-i to denote the irreps [Aat-i, AAr_ 2 ]. Note that 
the sum in Eq. (^^ over the irreps A'(r_^ is restricted to those states which belong both to 
the irrep Xn and X!j^. The kinematical “relative coordinate” rotation applied to the Jacobi 
coordinates 


^N-l = 


— 


2) rh-i + ta 



A-2 


A-2 

E o. 


i=l 


( 30 ) 


(29) 
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and given by 


= ~\l 

is realized in Eq. (|^) by the generalized, On, Wigner D fnnctions, where gA,A-i is the gronp 
element that corresponds to the rotation (|3lD , i.e. a rotation by an angle 7 = arcsin(Y^^^^^^) 
in the fjN,gN-i plane. 

Since the two-body operator in Eq. (p9D depends only on the coordinate ^n we shonld 
separate the hyperspherical fnnctions related to this coordinate in the bra and the ket states. 
So, the last step in the derivation consists of nsing the ohscfps introdnced in Eq. (pS]). Then, 
the last term on the rhs of Eq. is 


{Kn-iLn-iXn-iI3n-ii^n)KnLn\}KnLnXnI3n 

1 ijv- 1 /3^_ jv/k 

\Kn-iLn-i\n-i(3^-i,(.'n)K'nL'n\]K'^L'^X^(3^'^ 


X 


X < {Kn-iLn-i ■ ■ ■',^n)KnLnMn\0^{'\/2^n)\{Kn-iLn-i ■ ■ ■',^'n)^'nL'n^'n > ■ ( 32 ) 


This matrix element can be calcnlated for any given two-body operator. For scalar operator 
0^{rA,A-i), one obtains 


OKNK'^-,KN-itN^P) - 

< {Kn-iLn-1 .. .■,EN)KNLNMN\0^{V2iN)\{KN-iLN-i ... >= 

5,,,e,MN{KN, ^NKN-l)MN{K'J,■, ^NKN-l) [" d9sin^^^(9) 

J 0 

xPn {cos29)F^, ^ ^ {cos29)0 {v2psm9), 


(33) 


with n = = 


and where 


Mn{K_ 


N] rv 


Kn-i) = 


{2Kn + 3N- 2)n!r(n + Kn-i + In + 


r(ri + En + |)r(Ti + Kn-i + 


3V-3^ 


(34) 


is a normalization constant, and are the Jacobi polynomials. 

The spin-isospin part of the two-body operator is calcnlated in a similar way, hrst we nse 
the last two generations of the spin-isospin cfps ( Eq. (P2|)) to get the explicit dependence of 
the state > on the spin and the isospin of the last two-particles. 


< SAS^ATAT^ATAYA-iaf\0^^{sA,SA 
^Ya-2,Y'a.2 Y. Y. 


_1, tA, tA-l)\S'j,S'^AT'/r'^AT'AY'A.i aT >= 

[{SA-l,SA)SA{TA-i,tA)TATA-iaf_^}SATATACif] 


SA-iTA-iaPL 


, T' 

A—1 ^ A—1 A—1 *^A—1 


X [(s:,_i; S^)s;4(n_i; tA)T'A'A-l<^A-, 


\)S'aTA> 


/5Ti 

A J 
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X 


^ [{SA-2]SA)SA-l{TA-2]tA)TA-lTA-2af_^\}SA-lTA-lTA-iaf_^] 
Sa-2Ta-20l^j^_2 

X < SaTaSa-iTa-iSa-2Ta-2 ■ ■ ■ \0^^(SA, SA-1, Ia, tA-2)\S'AT'AS'A-2l"A-2S'A-2T'A-2 ' ' 


(35) 

. > . 


Then the matrix element on the rhs of (^) is calculated using the 6j symbols. 

Special attention should be payed to the matrix elements of the generalized angular 
momentum operator K as it contains the hyperangular part of the Laplace operator. Using 
the dehnition in Eq. (|T^ we see that 

Tk'u'p',Kup=< K'i''\xp'{zij)*kxp{zij)\Ki' > . (36) 


This matrix element can be easily calculated once we know the action of K on Xp(%)- 
Since K is invariant under kinematical rotations we can always write it using a set of Jacobi 
coordinates such that — ^j)) which leads to the following expression 


< K'v'\xp'{zijykxp{zij)\Kv >= ^ 

< K'v'\xA^A’Ki.K + 3Af - 2)xMi) “ X,'(2«)'(l - 


-2xp'{zij) (1 z^j) 


2\^Xp{^ij') ^ 


dzij dzij 


dz; 
\Ku > . 




(37) 


This is just the matrix element of a sum of two body operators which can be calculated as 
explained above. 

roj 

A bit more problematic are the three- and four-body integrals and 

p'Kup- order to calculate the three- (four-) body term we must first use three 
(four) generations of the ocfps in order to express the matrix elements in terms of the or¬ 
thogonal symmetry adapted hyperspherical functions. Then we can use the proper rotations 
and reduce the integrals to six-dimensional integrals that depend only on the last two Jacobi 
coordinates. The explicit dependence of the HH functions on the Jacobi coordinates is then 
revealed using the last two generations of the ohscfps, Eq. (|^. For L = 0 states, these 
integrals can further be reduced to three-dimensional integrations. 


V. DISCUSSION OF THE RESULTS 


The present method becomes more and more complex as the number of fermions in¬ 
creases. Therefore in this work we consider only central NN potentials. We present results 
for 6, 8 and 12 nucleons interacting via the Volkov |]^ (VV), the Afnan-Tang (S3), 
the modihed S3 potential (MS3), the Brink-Boeker [Q (Bl) and the Malfliet-Tjon 


2^ potentials (MT-I/III and MT-V). The accuracy and convergence of the method have 


been investigated for 6 and 8 nucleons. 

The first step in our numerical study was to determine the hyperradial working point, 
i.e. the parameters a and b of the Laguerre polynomials in Eq. (^OD , the number of radial 
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grid points and their location. After some preliminary tests we decided to use 15 hyperradial 
grid points and 15 hyperradial functions with the parameters b = 0.25 fm and a = 12. The 
grid points where chosen as the abscissas for a Gauss-Laguerre integration. 

The next step in setting our working point was to check the convergence of the two- 
body correlation terms. As an example we studied the ®Be system with the hard core 
MT-V potential, using a single HH function, with K = 4 and expanding the correlation 
function into Jacobi polynomials. As can be seen from Fig. the binding energy converges 
very fast with increasing number of polynomials . In fact, the energy difference between 
the calculations with 7 and 8 correlations terms is about 0.02 MeV. As a consequence we 
used 8 terms in the expansion of the correlation function. 

The effect of introducing the hypercentral potential I4c(p) is presented in Table |for ®Li 
with the VV potential. As can be seen the hypercentral potential accelerates the convergence 
of the HH expansion. With the lowest K^ax of 2 one already obtains a larger binding energy 
than with K^ax = 6 without hypercentral potential. Thus the introduction of Vhc is certainly 
advantageous. 

Our numerical results with the parameters described above are presented in Table for 
the 6 nucleon system and in Table 0 for the 8 nucleon system. 

For ®Li the calculations include only one line irreps of the kinematical group and the 
irreps [42] and [33] of the permutation group. We compare our values with recent accurate 
variational results [0 available for some of these central potentials. One can notice that 
with the VV potential one reaches convergence faster than with other potentials. The result 
for the binding energy starts oscillating around the asymptotic value. It also compares nicely 
to the variational result of Ref. [^. The other potentials show a tendency to convergence 
even if K^ax is not large enough to reach it. The MT potentials seem to lead to a more 
rapid convergence than S3 and Bl. The differences may be due both to the fact that the 
convergent value has not yet been reached and to the missing irreps of the permutation 
group. The two line irreps of the orthogonal group are of little importance. This has been 
checked for the 8 particle case, where they give rather small contributions. 

The calculations for ®Be include only the irrep (400) of the kinematical group Oj and 
irreps with at most 3 rows of the permutation group. Here the comparison is made with the 
TICI results [Q. In all cases our results for the binding energy are somewhat larger. They 
show characteristics similar to the six-body case. Again one sees that the VV potential 
result presents small oscillations around the convergent value and that the rather hard core 
MT potentials seem to give values closer to convergence than Bl or S3. From the comparison 
between the TICI results with and without state dependent correlations one can infer that 
for the MT potentials state independent correlations already give rather satisfying results, 
while state dependent correlations are more effective for Bl and S3. Since our correlations 
are state independent one could expect such a different convergence behavior as found in 
Table ^ 

The results for are presented in Table |I^. Here the calculations include a single HH 
state, Kmax = 8, with the kinematical group On irrep (4,4,0) and the 5i2 irrep [444]. The 
comparison is made with TICI [§, IDEA 0] , and the variational result of Ref. |p4||. Again 


from the comparison with the TICI results one can see that S3 and Bl results are farer from 
convergence than the MT values and that state dependent correlations play a similar role 
as discussed for the ®Be case. 
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The IDEA results for 12 particles deviate somewhat from ours, although they are all 
obtained using a single HH state. One possible explanation of the difference can be attributed 
to the choice of the HH state which is not unique. On the other hand we have obtained 
excellent agreement with the IDEA results for 16 bosons. 


VI. CONCLUSIONS 

In this work we have formulated a general method to calculate the wave functions of 
light systems up to considerably large number of particle. This method combines the main 
ideas of the HH expansion, the pair correlation ansatz and the Faddeev approach. The 
actual application of it is made possible by the use of a very efficient recursive algorithm 
to construct the antisymmetric A-particle state containing hyperspherical harmonics. We 
have applied the method to calculate the binding energies of 6, 8 and 12 nucleon systems 
with central local potentials. The results we have obtained are very encouraging. For some 
potentials (VV, 6 particles and MTV, 8 particles) we have reached the convergence region 
with Kmax = 8 which is the maximum value allowed by our present computer facilities 
(workstations). The ®Li result for the VV potential is slightly higher than the SVM result. 
Even if in other cases we have not yet reached the convergence in the HH expansion our 
results are close to the TICI results. For the eight-body case they are higher for all the 
potentials where results were available for a comparison. 

The method presented here for the solution of the few-body Schrodinger equation can 
be easily extended to solve the Schrodinger-like equation with a source, necessary for the 
application of the Lorentz integral transform method. Work in this direction is in progress. 
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TABLES 


TABLE I. 

The effect of the hypercentral potential Vhc on the binding energy of six-nucleon system ^Li, 
{L,S)J'^ = (0,1)1"'’, interacting via the VVpotential. 


K 

^^max 

Vhc = 0 


2 

64.18 

66.10 

4 

64.81 

66.53 

6 

65.47 

66.63 


TABLE II. 

Binding energies of six-nucleon system ^Li, {L,S)J'^ = (0,1)1^? interacting via various NN po¬ 
tentials. Nhh represents the number of hyperspherical harmonic states. 


K 

^^max 

Nhh 

B1 

MT I/III 

MTV 

S3 

VV 

2 

1 

30.99 

30.15 

62.45 

62.76 

66.10 

4 

4 

37.82 

34.67 

63.27 

64.45 

66.53 

6 

12 

39.11 

35.43 

64.10 

66.49 

66.63 

8 

31 

39.61 

35.91 

64.55 

67.18 

66.57 

SVM 1 

1 

- 

- 

66.30 

70.65 

66.25 


TABLE III. 

Binding energies of eight-nucleon system ^Be, {L,S)J'^ = (0,0)0'’', interacting via various NN 
potentials. Also given are results from Ref. ^ with state independent TICIsi and state dependent 
TICI SD correlations. 


K 

^^max 

Nhh 

B1 

MT-I/III 

MTV 

S3 

VV 

4 

1 

56.71 

52.82 

134.29 

31.19 

147.42 

6 

4 

65.39 

59.31 

137.72 

38.08 

148.70 

8 

15 

70.03 

60.64 

137.80 

42.11 

148.49 

TICI 5 / 


49.18 

46.67 

129.25 

26.26 


TICIsd 


61.30 

52.67 

130.23 

37.30 
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TABLE IV. 


Binding energy of twelve nucleon system 
potentials. Also given are results from Refs, 
results as in Table Wj . 


C, {L,S)J'^ = (0, interacting via various NN 
^ (TICI), ^ (IDEA) and ^ (VMC). The TICI 



B1 

MTM/III 

MTV 

S3 

VV 

This work 

96.64 

108.34 

437.25 

52.42 

494.00 

TICI 

103.93 

109.04 

429.44 

62.99 


IDEA 

80.1 



44.4 


VMC 

82.9 ± 0.2 
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